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Abstract 

This is a technical note where we solve the additive eigenvalue problem associated to a dy- 
namics of a 2D-traffic system. The traffic modeling is not explained here. It is available in [2]. It 
consists of a microscopic road traffic model of two circular roads crossing on one junction man- 
aged with the priority-to-the-right rule. It is based on Petri nets and minplus algebra. One of 
our objectives in [2] was to derive the fundamental diagram of 2D-traffic, which is the relation 
between the density and the flow of vehicles. The dynamics of this system, derived from a Petri 
net design, is non monotone and additively homogeneous of degree I. In this note, we solve the 
additive eigenvalue problem associated to this dynamics. 

1 Introduction 

In this note we solve the additive eigenvalue problem (or the time-independent system) associated 
to the dynamics of a basic 2D-traffic model considered in [2]. It is a system of two circular roads 
crossing on one junction managed with the priority-to-the-right rule. The model is based on Petri 
nets and on minplus algebra pQ. It is an extension to an existing ID-traffic model [31 [4], which gives 
the average speed of vehicles on one circular road as an eigenvalue of a minplus matrix, and thus 
allows the derivation of the fundamental diagram of ID-traffic (the relation between the density and 
the flow of vehicles on the road). 

We give a solution of the eigenvalue problem. We show that the eigenvalue is not necessarily 
unique, but is given in terms of two main quantities which are interpreted in terms of traffic as the 
density d of vehicles in the system, and a parameter r giving the ratio of the non priority road size 
with respect to the size of the whole system. We give a condition on the parameter r such that the 
eigenvalue is unique and positive for non- high densities. In this case the eigenvalue problem can be 
written as a dynamic programming equation of a stochastic optimal control problem. 

We use the minplus algebra notations, mainly for reason of compactness but also to use some 
classical results of this algebra p~j. In addition, the following notations are also used: afb denotes 
a — b, y/a denotes a/2, and b a denotes ab. 
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The traffic dynamics is the following (see [2]): 
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For example, in the usual algebra, the equation ([5]) is written: 



when the equation ([4| is written: 
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This system of equations is implicit but it is triangular, so its trajectory is unique. 

We denote by d the following quantity (which is interpreted in terms of traffic as the density of 
vehicles in the system): 
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2 Solving the additive eigenvalue problem 

The additive eigenvalue problem corresponding to the dynamics (DS) is: 

Xx l = Oj-iXj-i © OiXi+i, 

i £ {2, . . . , n — 1, n + 2, . . . , n + m — 1} , 
\x n = a ri xix n+ if{Xx n+m ) © a n -ix n -i , 

\EV) . \Xn-\-m — G n _)- m Xi fx n © Q, n + m —iX n -\- m — 1 , 

with the constraints ^ and the notation ([7]). Our aim in this note is to solve the system (EV) 
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Theorem 1. Solving the system (EV) is equivalent to solving the following simplified system (SS): 

Xi = (a,i-if\)xi-i © {a,if\)x i+ x, (13) 

i £ {2, . . . , ri — 1, n + 2, . . . , n + m — 1} , 
x„ = (a„/A 2 )xix„ +1 /a; n+m (6„/A™^ 1 )xi , (14) 
(SS) : /A)zi © (b m fX n ~ 1 )xn+i , (15) 

aci = (a n f\) v / x n x n+m © (6 n /A" _1 )a;„ , (16) 

(17) 



where b n = 0" = / a. t , 6„ = ®" =1 a;, & m = (g>™^„™ / a, and 6 m = (g)™ 4 



<l,. 



Proof. We proceed in two steps: 



• First we show that if (A, x) is a solution of the system (EV) then A < 1/4. Indeed, from the 
equations ©, (fTTj) and (TT!?j) . we have: 

Ax„ = a„a;ix Tl _)_i/(Ax ri 4_ TO ) © a„_ix„_i, 

Axi $n ^/ XnXn+m : 
AXn+1 fln+m^ X n X n -^ m . 

Then by multiplying (standard adding) the terms of these inequalities, we obtain A 4 < 1, since 

Q'Ti+mQ'nQn+m 1- 

• We see that if n = m = 2 the systems (EV) and (SS) coincide. For n and m fixed, we 
denote by EV(n, m) and SS(n, m) the corresponding systems. By induction on n and m, we 
suppose that EV(n,m) SS(n,m), and we show that EV(n + 1, m) 4=> iSS^n + l,m) and 
EV(n, to + 1) S^n, m + 1). 

— To show that EV (n + 1, to) SS(n + 1, m), we eliminate the variable .t„ in EV(n + 1, m) 
which gives a system EV(n,m), then we use the induction assumption. 
Indeed the problem EV(n + 1, to) is written as follows: 

XXi = CLi-xXi-X © OiXi+i, (18) 

i G {2, . . . , n, n + 3, . . . , n + to} , 
Ax„+i = a„ + iXia; Il+ 2/(A2:„ + i +m ) © a n x n , (19) 

Ax n -|_l-|- m Q Jn J r \J rm X\X n J r 2fx n J r \ © Q, n + rn X n -\- m , (20) 
Axi = On+i^/Xn+Un+i+m © ai2;2 , (21) 
Ax n +2 = fl|i+l+mv'5i+Si+l+m © 5n+2^n+3 7 (22) 

Using the expression of x n in (|18p . we replace it in the expression of in (|19p . We 

obtain: 



Ax„+i = a„ + iXia;„ + 2/(Aa;„ + i +m ) © a n [(a„_i/A):r„_i © (a„/A)x 



n+ljj 
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which gives: 

Ax„+i = a n+ ixix n+2 f (\x n+ i +m ) © (a n a n -ifX)x n -i, (23) 
because Xx n+ i < (a n a n fX)x n+i since A<l/4<l/2 and a n a n = 1- 

Using the expression of x n in (JTHJ) , we replace it in the expression of x n -\ in (JT5J) also. We 
obtain: 

Az„_i = a„_ 2 a;ri-2 © 5„_i[(a„_i/A)a; Tl _ 1 © (a„/A)x n+ i], 

which gives: 

Ax„_i = a n - 2 x n -2 © (a n _ia„/A)a; n+ i, (24) 

because Acc„_i < (a„_ia„_i/A)a; rl _i, since A < 1/4 < 1/2 and a n _ia„_i = 1. 
Let us denote by X\ and X 2 the following sets of indexes: 

Xi = {l<i<n + 1 + m, i =/= 1, n + 1, n + 2, n + 1 + to}, 

X2 = {1 < z < ?i + 1 + to, i ^ 1, 71 — 1, 71, » + 1, 71 + 2, ri + 1 + to}. 



Thus we can conclude the following equivalence: 

& TOieij, mm, gb><=«]. (25) 

The equations |JTH])iex 2 combined with the equations pip. (|23|) . (|20|) . (f2~Tj) and ([22]) form 
the following EV(n, to) system of variables Xi, l<i<n + l + m and i =/= n: 

Xxi = a^x^i © a,iX i+1 , (26) 

i G {2, . . . , n — 2,n + 3, . . . , n + to} , 
Ax„_i = a n ^ 2 x„-2 © (a„_ia„/A)x„ + i , (27) 
Ax„ + i = a„+ia;ia;„ + 2/(Aa;„ + i +m ) © (a n a n -ifX)x n -i , (28) 

(29) 

Axi = On+i^/Xn+iIn+i+m © 5iX 2 , (30) 

(31) 

By using the induction assumption, this later system is equivalent to the following system: 

Xi = (a i -if\)x i - 1 © (a l fX)x i+u (32) 
i£ {2, . . . ,n — 1, 7i + 3, . . . , 7i + to} , 

/A ) © (&„+i/A n )a:i , (33) 

= (an+i+mfX)xix n+ 2fx n+ i © (b m fX m 1 )x n+2 , (34) 

El = (a n+ i/A)^x„+ia;„ + i +m © {b n +ifX n )x n +i , (35) 

*£n+2 (^n+l+m 

fX)y/x n+1 x n+1+m © (b m fX m ) (36) 



By adding the equation (|18p ,-n to the later system, we obtain SS(n + 1, to). 
— We show with the same manipulations that EV(n, to + 1) SS(n 7 m + 1) ■ 
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Lemma 1. (Baccelli et al. fJjj) Given A a (m x m) minplus matrix, if the weights of all the circuits of 
the graph Q(A) associated to A are positive, then the equation x = A®x@b admits a unique solution 
x = A* ®b, where 



A*=0A" = A r ' 



n=0 



n=0 



Corollary 1. Solving the system (EV) is equivalent to solve the following system of four variables: 

x„ = (a n f\ 2 )x 1 x n+1 fx n+m © (6„/A" _1 )a;i , 
(S) : x n+m = (a n+m f\)xix n +ifx n © (b m f\ m ~ 1 )x n+ i , 

x\ = (a„fX)^x n x n+m © (b n fX n ~ 1 )x n , 



(37) 
(38) 
(39) 
(40) 



Proof. Taking into account the equivalence (EV) (SS), suppose that A, X\, x n , x n+ i, x n+m are 
known. To determine the other variables i.e. X2, • ■ ■ , x n —i, x n +\, . . . , x„ +m _i we have to solve the 
system of equation l|8l). This system is written: 



(41) 



where 



x — (x% , . . . , x n —\ j x n .^_i , . . . , x n _j_ m _i ) , 



with 



A= (e/A) 

e 0,2 £ 
a 2 £ a 3 

e a 3 e 
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£ £ 
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. . . , e 



^2 = 



and 

b = t [aiXi , e , 

Using Lemma [T] all the circuits of the graph associated to the matrix A have the average weight 1/A 2 
which is positive because A < 1/4 < 1/2. Thus the matrix A* exists and the solution of the system 
([4T|) is given by: x = A* <g> b ■ 
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3 Solving the system (S): 



Let us use the notations: 

r = nj (n + to — 1) , 

p = l/(n + m — 1) = r/n . 
Theorem 2. There exists a solution (A, a;) of S such that A satisfies: 

= max jmin jd - (1 + p)A, i - A, r - d - (2r - 1 + p) A j , -A 

Remark 1. Before we give the proof of Theorem [21 let us explain it. Using the notations: 

di = (n + m)/[4(n + m - 1)] = (1 + p)(l/4) , 

d 2 = (3n + m- 2)/[4(n + m- 1)] = (2r + 1 -p)/4 , 
the result can be explained as follows (see Figure [1]): 

• If < d < d\ then (S) admits a solution (A, x) such that: X = d/(l + p) , 

• If di < d < efe the (S) admits a solution (X,x) such that: A = 1/4 . 

• If c?2 < ^ < r or r < d < d 2 which cases correspond respectively to r > 1/2 or r < 1/2 then (S) 
admits a solution (A, x) such that: A = (r — d)/(2r — 1 + p) , 

• If r < d < 1 then (S") admits a solution (A, x) such that: A = . 
Proof. 

• If < d < di , then a solution (A, x) is given by: 

d n + to — 1 
A = = a, 



1 + P 



n + m 



x n 




" bn/A"" 1 






A" +1 /<4/&„ 


si 




e 






@ / n-\-m/@"n 



which is a solution of: 



Indeed: 
[(a„/A 2 ) 



= (Kf A™ 

^n+m (^m/A )-^n+l; 
3?n+l : — (^n+m/A)-^/ Xn%n-\-r 



1/A 4 , because a n = l/(a n a n+m ), 



> e, because A = d < ^T 1 " 1 1/4 < 1/4. 



G 



[(a n +m/A)xix„ + i/x n ] /x„ +m = 1/A 3 , because a n+m = l/(a„a„ +m ), 

> e, because A < 1/4 < 1/3. 
[(&„/A n " 1 )x„] f Xl = (1/A 2 )"- 1 , because b n b n = l"-\ 

> e, because A < 1/4 < 1/2. 
[(&m/A m_1 )x„ +m ] /x„ +1 = l m - 1 /(6 Il 6 m a Il a„ +rn )A"-™+ 2 , because 6 m = l™" 1 /^, 

= (1/A 2 )" 1 - 1 , because b n b m a n a n+m = d n+m - x = \ n + m , 
> e, because A < 1/4 < 1/2. 

If g?i < d < c?2 then a solution (A, x) is given by: 

A = 1/4 







A™- 3 a„/6 m 








Xl 




e 


.Xn+1. 







which is a solution of: 

Xn (^n/A )xiX n +i/x n _|_ m , 
Xn+m — (^m/A )x n +i, 
Xl {p>nf A)-^/x n x n ^ m , 
Xn+1 — (^n+7n / A) tj%n ^n-\-m • 

Indeed: 

[(&n/A" _1 )xi] /x„ = (a n a„ +m 6 n 6 m )/l/A" +m - 4 , because a„ = l/(a n a n+m ), 

= rf"+ m - 1 /l/(l/4)" +m - 4 , 

because a n a n+m b n b m = d" +m " 1 , and A = 1/4, 
> B±m _ x _ n+rpi = ^ because d > dl . 



[(a„ +m /A)xix„ + i/x„] /x„ +m = A = 1/4 > e. 
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[(pnfX n ~ 1 ) x n\ fxi = l n f{a n a n+m b n b m )fX m - n -\ 

because b„ = l n ~ 1 fb n , and a n = \f(a n a n+m ), 

= l«/d«+™-i/A"- m + 2 , because a n a n+m b n b m = d n+m ~ l , 

> n - 3n +™- 2 - "-^+ 2 = e , because d < d 2 . 
[{l m f\ m - x )x n+m ] fx n+1 = (l/A 2 )™" 1 , because b m b m = I™- 1 , 



> e, because A = 1/4. 



• If d 2 < d < r or r < d < d 2 , then 

A 



r — 


d 


n 


n 


+ m 


- 1 


2r - 1 


+ P 


n — m - 


- 2 n 


— m 


+ 2 






A™ 


-'fbn 










b m a n+mt 




-n+l 




Xi 






e 




1 


_X n +l_ 




_ b m dn +m f 


b n f\ m - 


n+2 





which is a solution of: 

x\ = (6 n /A™ _1 )o;„ , 
Indeed, A given by (|42|) satisfies < A < 1/4 because: 



(42) 



rt-f-m— 1 



> 0, and we can check 



1. If da < d < r, which corresponds to n - m + 2 > 0, then n 
that: d 2 <d<r => < A < 1/4. 

2. If r < d < d 2 , which corresponds to n — m + 2 < 0, then < 0, and we can check 



that: r < d < d 2 => < A < 1/4. 



Then: 



i n-l 



[(&n/A" _1 )a;i] fx n = [I/A 2 ]' 

> e, because A < 1/4 < 1/2. 



i n—1 
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[(a n+m f^)XlX n+1 fx n ] fXn+r 



n \m—n—l . 



1"A 



a n fa n+ ifb m fb n , 



[(a n f\) y /x n x n+m ] fxi 



[(b m fX m l )x n+m ] fx n+1 



because a n+m = lfa„fa n+m and b n = 1" 1 fb n , 



= ( ' 1 « A r l -m+2 / ( d n+m-l- )A ^ 

because a n a n+m b m b n = 1 

= A > e , because A"- m + 2 = l n fd n+m - 1 . 
= a„a„ +m 6„6 m /l n_1 A n " m ~ 2 , because b n = \ n ^ x fb ni 

= (d"+ m - 1 /l n A"-" l+2 )(l/A 4 ) , because a n a n+m b m b n = d n+m - 1 

= 1/A 4 because A""™ 4 " 2 = l"/^ 4 "" 1 " 1 , 

> e, because A < 1/4. 

= i'«-i/A 2m - 2 , because b m b m = l™" 1 , 

= a/A 2 r-\ 

> e, because A < 1/4 < 1/2. 



If r < d < 1, then A = 







efb n 


Xn-\-m 




V l+1 ta 2 Jb n 


X\ 




e 


_X n +l _ 







which is a solution of: 



Indeed; 

[(bn?X n ~ 1 )x 1 ] fx r 



x n - {o, n f\ ^X±X n ^.\ f X n -\- m , 

x\ = (b n f\ n ~ l )x n , 

Xn-\-l {(Ln+mftyy/ X n X n -\- ra , 



l n 1 >e, because 6 n 6 n = l n , 



[(^m/^ )^n+l] f^n-\-m ^n^n+m / 1 j 

= rf»+ m -Vl™, because a n a n+m b n b m = d n+m ~\ 
> e, because d > r. 
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[{a n fX) y /x n x n+m ] fxi = 1 >e, 

[(b m f\ m ~ 1 )x n+m ] fx n+ i = l^" 1 - 1 f(a n a n+m b n b m ) , because b m = l m ^ 1 fb m , 

= i«+m-iy d n+m-i ) because a„a n+m b n b m = d n+m ~ l , 

> e , because d < 1 ■ 

Corollary 2. In £/*e case r > 1/2, a non negative eigenvalue A o/ (5) is given by: 

f . f 1 ,1 r-d 1 

A = max < mm < a , - , > , 

\ \ 1 + p ' 4 ' 2r - 1 + p J ' 

Proof, follows directly from Theorem [2] ■ 

Corollary 3. For large values of n and m such that n > m — 2 (which is the case r > 1/2), a non 
negative eigenvalue A of (S) is given by: 

A = max | min Id, —, 1 , 

\ \ ' 4 2r - 1 J 

Proof, follows directly from Corollary [3] ■ 

Remark 2. We can check that as soon as we assume m > 1, we get d\ < d 2 , so we have: < d\ < 
d 2 < 1- The position of r with respect to d\ and d 2 gives three cases and divides the interval [0, 1], in 
each case in four regions, where the eigenvalues A of (S) satisfy (see Figure [1]): 

A. d e [0, min(di, r)[ => A = -p- d, 

{or A = ^r 5 d_ i , 
or A = 0. 

C. d £ [di,min(d 2 ,r)[ => A = ±, 

A = -, 

D. d G [max(di, r), d 2 [ '(or A 



r— <i 



27'-! + ,!) ' ' 

or A = 0. 

E. de [d 2 ,max(d 2 ,r)[ => A = 2r r _~+ p , 

F. d g [max(d 2 ,r), 1] => A = 0. 

Lemma 2. If X > then the system (S) (solved on (A, x) j is equivalent to the following system (solved 
on (A, z) ): 

z„ = (d n ?b m fX 2 ) Zl © (b n fX n+m - 2 ) Zl , (43) 
(SZ) : z n+m = b m z n+1 , (44) 

zi = (a n f\)^/z n z n+m ® (6„/A"" m )z„. (45) 
z„ + i = (a„ +ln /A) v /z„2:„ +m © (& m /A 2m ~ 2 )z n+m . (46) 
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1/4- 



B D 



1/4- 



1 



A C 



1/4- 



CASE 1 



d, r d 2 

CASE 2 



1 o 



A 


C 


E 


F 


/////// 1 






d 



CASE 3 



Figure 1: The curve of A given in Theorem [2] depending on d. 



Proof. From the equations (|57|) and ([55)) we obtain: 

XnXn+m (^n+m/A)xiX n -|_i . 

Since A > we have: x n x n + m < (a n + m f\)xix n +i. Thus: 

x n+m = (b m f\ m - 1 )x n+1 . (47) 

By replacing x n+m in (|37[) . we obtain: 

x„ = [a„/6 m A m - 3 © 6,,/A"- 1 ] x\ . (48) 

The system (S*) is then equivalent to the system { ([4*51) . (|4"7) . (j3"5j) . (|4*U| }. On this later system we use 
the following changing variable: 

_ \2m-2 
Zn+m — Xn-\-mA , , 

Zri+1 = X„ + iA m , 

and we obtain the system (SZ) ■ 

Theorem 3. If r > 1/2 (that is n> m), and for the densities d satisfying < d < r, the system (S), 
and thus the eigenvalue problem (EV), admit a unique positive eigenvalue A given by: 

A = min{— i— d, \ , - r ~ d } >0. (49) 

This situation corresponds to the phases A, C and E of the case 3 on Figured of Remark® 
Proof. 

• Let A be positive. Lcmma[2]gives the equivalence of the systems (S) and (SZ). The later system 
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is the eigenvalue problem associated to the following dynamical system: 

z k n = (a n fb m )zt 2 © b n z k - [n+m - 2 \ (50) 

z n+m = b m Z n+ i, (51) 
z\ = a n \j ' Zn 1 z k + l n © bnZn ^" m \ (52) 

*n+l - »»+m^ ^n+m w "m z n+m ■ l 00 / 1 

If r > 1/2, that is n > m, then this dynamical system is implicit but triangular. Indeed, an 
iteration of the dynamics is to compute z k and z k +1 in parallel, then compute z k +m , and finally 
compute z\. So the system (SZ) can be interpreted as a dynamic programming equation of a 
stochastic optimal control problem, where A is the average optimal cost by unit of time. Since 
A is supposed to be positive, and from Corollary [21 we obtain ([49]) . 

• Let Ai > and A2 > be two positive eigenvalues of (S). Lemma [2] tells us that both of Ai 
and A2 are eigenvalues of the system (SZ). Since in the case when r > 2, the system (SZ) is 
a dynamic programming equation of a stochastic optimal control problem, which thus admits a 
unique eigenvalue, we conclude that Ai = A2 ■ 

4 Conclusion 

The results of this note give a solution to the additive eigenvalue problem associated to a dynamics of 
an elementary 2D-traffic system (two circular roads crossing at one junction, managed by the priority- 
to-the-right-rule) . The eigenvalue A, which is not necessarily unique, is given as a function of two 
main quantities which are interpreted in terms of traffic as the density d of vehicles in the system, and 
the ratio r between the non priority road size and the size of the whole system. Moreover, we showed 
that when r satisfies r > 1/2, that is when the size of the non priority road is bigger than the size of 
the priority road, the uniqueness of A, which is positive in this case, is proved for densities satisfying 
d < r. 
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